(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     33998,        944]
NotebookOptionsPosition[     32426,        888]
NotebookOutlinePosition[     32856,        905]
CellTagsIndexPosition[     32813,        902]
WindowFrame->Normal
ContainsDynamic->False*)

(* Beginning of Notebook Content *)
Notebook[{
Cell["\<\
First define the Antenna Pattern Functions. Note that A and B (as defined in \
JKS 98) are functions of sky - cordinates (\[Theta], \[Phi]). Also note that \
\[Psi] is the polarization angle.\
\>", "Text",
 CellChangeTimes->{{3.472931937116972*^9, 3.472931976053673*^9}, {
  3.472932046406058*^9, 3.47293249090398*^9}, {3.472932531127093*^9, 
  3.47293262819089*^9}, {3.472932682057527*^9, 3.472932705260825*^9}, {
  3.472932759786611*^9, 3.472932767263685*^9}, {3.472932805887068*^9, 
  3.472932923089784*^9}, {3.472933002871185*^9, 3.472933014102472*^9}, {
  3.47293305532964*^9, 3.472933076145941*^9}, {3.472933158541131*^9, 
  3.472933191015629*^9}, {3.472933297221234*^9, 3.472933362270807*^9}, {
  3.472933405647537*^9, 3.47293344020603*^9}, {3.472933630455057*^9, 
  3.472933632371549*^9}, {3.472933760549579*^9, 3.472933988294928*^9}, {
  3.472934022424597*^9, 3.472934068201747*^9}, {3.472934116665703*^9, 
  3.472934132496888*^9}, {3.472934166636935*^9, 3.47293428014505*^9}, {
  3.472934440035422*^9, 3.472934440330311*^9}},
 FontSize->18],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   SubscriptBox["F", "p"], " ", "=", " ", 
   RowBox[{
    RowBox[{"A", " ", 
     RowBox[{"Cos", "[", 
      RowBox[{"2", " ", "\[Psi]"}], "]"}]}], " ", "+", " ", 
    RowBox[{"B", " ", 
     RowBox[{"Sin", "[", 
      RowBox[{"2", " ", "\[Psi]"}], "]"}]}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   SubscriptBox["F", "c"], " ", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"-", "A"}], " ", 
     RowBox[{"Sin", "[", 
      RowBox[{"2", " ", "\[Psi]"}], "]"}]}], " ", "+", " ", 
    RowBox[{"B", " ", 
     RowBox[{"Cos", "[", 
      RowBox[{"2", " ", "\[Psi]"}], "]"}]}]}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.472931937116972*^9, 3.472931976053673*^9}, {
  3.472932046406058*^9, 3.47293249090398*^9}, {3.472932531127093*^9, 
  3.47293262819089*^9}, {3.472932682057527*^9, 3.472932705260825*^9}, {
  3.472932759786611*^9, 3.472932767263685*^9}, {3.472932805887068*^9, 
  3.472932923089784*^9}, {3.472933002871185*^9, 3.472933014102472*^9}, {
  3.47293305532964*^9, 3.472933076145941*^9}, {3.472933158541131*^9, 
  3.472933191015629*^9}, {3.472933297221234*^9, 3.472933362270807*^9}, {
  3.472933405647537*^9, 3.47293344020603*^9}, {3.472933630455057*^9, 
  3.472933632371549*^9}, {3.472933760549579*^9, 3.472933988294928*^9}, {
  3.472934022424597*^9, 3.472934068201747*^9}, {3.472934116665703*^9, 
  3.472934132496888*^9}, {3.472934166636935*^9, 3.47293430475754*^9}, {
  3.472934447295761*^9, 3.472934564616676*^9}, {3.473088734274578*^9, 
  3.473088858188752*^9}, {3.473088920979995*^9, 3.473088921648531*^9}},
 FontSize->16],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"Clear", " ", "[", "c", "]"}], "\[IndentingNewLine]", 
 RowBox[{"Simplify", "[", 
  RowBox[{"Sqrt", "[", 
   RowBox[{
    RowBox[{
     RowBox[{
      SubscriptBox["F", "p"], "^", "2"}], 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"1", "+", 
        RowBox[{"c", "^", "2"}]}], ")"}], "^", "2"}]}], " ", "+", " ", 
    RowBox[{
     RowBox[{
      SubscriptBox["F", "c"], "^", "2"}], " ", "4", " ", 
     RowBox[{"c", "^", "2"}]}]}], "]"}], "]"}]}], "Input",
 CellChangeTimes->{{3.474044190101524*^9, 3.474044321851211*^9}}],

Cell[BoxData[
 SqrtBox[
  RowBox[{
   RowBox[{
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{"1", "+", 
       SuperscriptBox["c", "2"]}], ")"}], "2"], " ", 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{
       RowBox[{"A", " ", 
        RowBox[{"Cos", "[", 
         RowBox[{"2", " ", "\[Psi]"}], "]"}]}], "+", 
       RowBox[{"B", " ", 
        RowBox[{"Sin", "[", 
         RowBox[{"2", " ", "\[Psi]"}], "]"}]}]}], ")"}], "2"]}], "+", 
   RowBox[{"4", " ", 
    SuperscriptBox["c", "2"], " ", 
    SubsuperscriptBox["F", "c", "2"]}]}]]], "Output",
 CellChangeTimes->{{3.474044217158686*^9, 3.474044281441236*^9}, 
   3.474044322372575*^9}]
}, Open  ]],

Cell["\<\
Now define the lowest order PN expansion terms (Sintes Vecchio). Note that we \
omit an explicit overall factor (2 M \[Eta] x)/r\
\>", "Text",
 CellChangeTimes->{{3.472931937116972*^9, 3.472931976053673*^9}, {
   3.472932046406058*^9, 3.47293249090398*^9}, {3.472932531127093*^9, 
   3.47293262819089*^9}, {3.472932682057527*^9, 3.472932705260825*^9}, {
   3.472932759786611*^9, 3.472932767263685*^9}, {3.472932805887068*^9, 
   3.472932923089784*^9}, {3.472933002871185*^9, 3.472933014102472*^9}, {
   3.47293305532964*^9, 3.472933076145941*^9}, {3.472933158541131*^9, 
   3.472933191015629*^9}, {3.472933297221234*^9, 3.472933362270807*^9}, {
   3.472933405647537*^9, 3.47293344020603*^9}, {3.472933630455057*^9, 
   3.472933632371549*^9}, {3.472933760549579*^9, 3.472933988294928*^9}, {
   3.472934022424597*^9, 3.472934068201747*^9}, {3.472934116665703*^9, 
   3.472934132496888*^9}, {3.472934166636935*^9, 3.47293430475754*^9}, {
   3.472934447295761*^9, 3.472934564616676*^9}, {3.473110580023403*^9, 
   3.473110584159238*^9}, 3.473132113034629*^9},
 FontSize->18],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   SubscriptBox["h", "p"], "=", 
   RowBox[{"TrigExpand", "[", 
    RowBox[{
     RowBox[{"-", 
      RowBox[{"(", 
       RowBox[{"1", " ", "+", " ", 
        RowBox[{"c", "^", "2"}]}], ")"}]}], " ", 
     RowBox[{"Cos", "[", 
      RowBox[{"\[CapitalPhi]", " ", "+", " ", 
       SubscriptBox["\[CapitalPhi]", "0"]}], "]"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   SubscriptBox["h", "c"], "=", 
   RowBox[{"TrigExpand", "[", 
    RowBox[{
     RowBox[{"-", "2"}], "c", " ", 
     RowBox[{"Sin", "[", 
      RowBox[{"\[CapitalPhi]", "+", " ", 
       SubscriptBox["\[CapitalPhi]", "0"]}], "]"}]}], "]"}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.472931937116972*^9, 3.472931976053673*^9}, {
   3.472932046406058*^9, 3.47293249090398*^9}, {3.472932531127093*^9, 
   3.47293262819089*^9}, {3.472932682057527*^9, 3.472932705260825*^9}, {
   3.472932759786611*^9, 3.472932767263685*^9}, {3.472932805887068*^9, 
   3.472932923089784*^9}, {3.472933002871185*^9, 3.472933014102472*^9}, {
   3.47293305532964*^9, 3.472933076145941*^9}, {3.472933158541131*^9, 
   3.472933191015629*^9}, {3.472933297221234*^9, 3.472933362270807*^9}, {
   3.472933405647537*^9, 3.47293344020603*^9}, {3.472933630455057*^9, 
   3.472933632371549*^9}, {3.472933760549579*^9, 3.472933988294928*^9}, {
   3.472934022424597*^9, 3.472934068201747*^9}, {3.472934116665703*^9, 
   3.472934132496888*^9}, {3.472934166636935*^9, 3.47293430475754*^9}, {
   3.472934447295761*^9, 3.472934588054855*^9}, {3.472934673065534*^9, 
   3.472934718164868*^9}, 3.473087882241709*^9, {3.473088270158252*^9, 
   3.473088277129054*^9}, {3.473088877333138*^9, 3.473088926159554*^9}, {
   3.473110593967132*^9, 3.473110624972783*^9}, {3.473131991770088*^9, 
   3.473131995799583*^9}, {3.473132046382698*^9, 3.473132048886281*^9}},
 FontSize->16],

Cell["\<\
Now expand the full waveform arriving at the detector, and collect terms as \
appropriate (See Ajith and Sukanta' s paper)\
\>", "Text",
 CellChangeTimes->{{3.472931937116972*^9, 3.472931976053673*^9}, {
  3.472932046406058*^9, 3.47293249090398*^9}, {3.472932531127093*^9, 
  3.47293262819089*^9}, {3.472932682057527*^9, 3.472932705260825*^9}, {
  3.472932759786611*^9, 3.472932767263685*^9}, {3.472932805887068*^9, 
  3.472932923089784*^9}, {3.472933002871185*^9, 3.472933014102472*^9}, {
  3.47293305532964*^9, 3.472933076145941*^9}, {3.472933158541131*^9, 
  3.472933191015629*^9}, {3.472933297221234*^9, 3.472933362270807*^9}, {
  3.472933405647537*^9, 3.47293344020603*^9}, {3.472933630455057*^9, 
  3.472933632371549*^9}, {3.472933760549579*^9, 3.472933988294928*^9}, {
  3.472934022424597*^9, 3.472934068201747*^9}, {3.472934116665703*^9, 
  3.472934132496888*^9}, {3.472934166636935*^9, 3.47293430475754*^9}, {
  3.472934447295761*^9, 3.472934588054855*^9}, {3.472934673065534*^9, 
  3.472934718164868*^9}},
 FontSize->18],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Collect", "[", 
  RowBox[{
   RowBox[{"TrigExpand", "[", " ", 
    RowBox[{
     RowBox[{
      SubscriptBox["F", "p"], 
      SubscriptBox["h", "p"]}], " ", "+", " ", 
     RowBox[{
      SubscriptBox["F", "c"], 
      SubscriptBox["h", "c"]}]}], "]"}], ",", "  ", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"A", " ", 
      RowBox[{"Cos", "[", "\[CapitalPhi]", "]"}]}], ",", " ", 
     RowBox[{"A", " ", 
      RowBox[{"Sin", "[", "\[Phi]", "]"}]}], ",", " ", 
     RowBox[{"B", " ", 
      RowBox[{"Cos", "[", "\[CapitalPhi]", "]"}]}], ",", " ", 
     RowBox[{"B", " ", 
      RowBox[{"Sin", "[", "\[CapitalPhi]", "]"}]}]}], "}"}], ",", " ", 
   "FullSimplify"}], "]"}]], "Input",
 CellChangeTimes->{{3.472931937116972*^9, 3.472931976053673*^9}, {
  3.472932046406058*^9, 3.47293249090398*^9}, {3.472932531127093*^9, 
  3.47293262819089*^9}, {3.472932682057527*^9, 3.472932705260825*^9}, {
  3.472932759786611*^9, 3.472932767263685*^9}, {3.472932805887068*^9, 
  3.472932923089784*^9}, {3.472933002871185*^9, 3.472933014102472*^9}, {
  3.47293305532964*^9, 3.472933076145941*^9}, {3.472933158541131*^9, 
  3.472933191015629*^9}, {3.472933297221234*^9, 3.472933362270807*^9}, {
  3.472933405647537*^9, 3.47293344020603*^9}, {3.472933630455057*^9, 
  3.472933632371549*^9}, {3.472933760549579*^9, 3.472933988294928*^9}, {
  3.472934022424597*^9, 3.472934068201747*^9}, {3.472934116665703*^9, 
  3.472934132496888*^9}, {3.472934166636935*^9, 3.47293430475754*^9}, {
  3.472934447295761*^9, 3.472934588054855*^9}, {3.472934673065534*^9, 
  3.472934718164868*^9}, {3.473087866535627*^9, 3.473087867084682*^9}, {
  3.473087917089607*^9, 3.473087917792158*^9}, {3.473110628851353*^9, 
  3.473110642561312*^9}},
 FontSize->16],

Cell[BoxData[
 RowBox[{
  RowBox[{"B", " ", 
   RowBox[{"Cos", "[", "\[CapitalPhi]", "]"}], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{
      RowBox[{"-", 
       RowBox[{"(", 
        RowBox[{"1", "+", 
         SuperscriptBox["c", "2"]}], ")"}]}], " ", 
      RowBox[{"Cos", "[", 
       SubscriptBox["\[CapitalPhi]", "0"], "]"}], " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"2", " ", "\[Psi]"}], "]"}]}], "-", 
     RowBox[{"2", " ", "c", " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "\[Psi]"}], "]"}], " ", 
      RowBox[{"Sin", "[", 
       SubscriptBox["\[CapitalPhi]", "0"], "]"}]}]}], ")"}]}], "+", 
  RowBox[{"A", " ", 
   RowBox[{"Sin", "[", "\[CapitalPhi]", "]"}], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"2", " ", "c", " ", 
      RowBox[{"Cos", "[", 
       SubscriptBox["\[CapitalPhi]", "0"], "]"}], " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"2", " ", "\[Psi]"}], "]"}]}], "+", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"1", "+", 
        SuperscriptBox["c", "2"]}], ")"}], " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "\[Psi]"}], "]"}], " ", 
      RowBox[{"Sin", "[", 
       SubscriptBox["\[CapitalPhi]", "0"], "]"}]}]}], ")"}]}], "+", 
  RowBox[{"A", " ", 
   RowBox[{"Cos", "[", "\[CapitalPhi]", "]"}], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{
      RowBox[{"-", 
       RowBox[{"(", 
        RowBox[{"1", "+", 
         SuperscriptBox["c", "2"]}], ")"}]}], " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "\[Psi]"}], "]"}], " ", 
      RowBox[{"Cos", "[", 
       SubscriptBox["\[CapitalPhi]", "0"], "]"}]}], "+", 
     RowBox[{"2", " ", "c", " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"2", " ", "\[Psi]"}], "]"}], " ", 
      RowBox[{"Sin", "[", 
       SubscriptBox["\[CapitalPhi]", "0"], "]"}]}]}], ")"}]}], "+", 
  RowBox[{"B", " ", 
   RowBox[{"Sin", "[", "\[CapitalPhi]", "]"}], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "2"}], " ", "c", " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "\[Psi]"}], "]"}], " ", 
      RowBox[{"Cos", "[", 
       SubscriptBox["\[CapitalPhi]", "0"], "]"}]}], "+", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"1", "+", 
        SuperscriptBox["c", "2"]}], ")"}], " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"2", " ", "\[Psi]"}], "]"}], " ", 
      RowBox[{"Sin", "[", 
       SubscriptBox["\[CapitalPhi]", "0"], "]"}]}]}], ")"}]}]}]], "Output",
 CellChangeTimes->{
  3.47308818913146*^9, 3.473088932050489*^9, {3.473110535024919*^9, 
   3.473110555370811*^9}, 3.473110604158832*^9, {3.473110643631525*^9, 
   3.473110661041683*^9}, 3.473130462364493*^9, 3.4731320010133266`*^9, 
   3.473132444888644*^9},
 FontSize->16]
}, Open  ]],

Cell[BoxData[
 RowBox[{" ", 
  RowBox[{
   RowBox[{
    RowBox[{
    "We", " ", "now", " ", "have", " ", "the", " ", "explicit", " ", 
     "expressions", " ", "of", " ", 
     SubscriptBox["h", "i"], 
     RowBox[{"(", "t", ")"}], " ", "and", " ", 
     SubscriptBox["a", "i"], " ", "where", " ", "the", " ", "full", " ", 
     "waveform", " ", "is", " ", "given", " ", "by", " ", "h", 
     RowBox[{"(", "t", ")"}]}], " ", "=", " ", 
    RowBox[{
     RowBox[{
      UnderoverscriptBox["\[Sum]", 
       RowBox[{"i", "=", "4"}], "4"], 
      RowBox[{
       SubscriptBox["a", "i"], 
       SubscriptBox["h", "i"], 
       RowBox[{
        RowBox[{"(", "t", ")"}], ".", " ", "Note"}], " ", "x"}]}], " ", "=", 
     " ", 
     RowBox[{
      SuperscriptBox[
       RowBox[{"(", 
        RowBox[{"\[Pi]MF", 
         RowBox[{"(", "t", ")"}]}], ")"}], 
       RowBox[{"2", "/", "3"}]], "has", " ", "been", " ", "absorbed", " ", 
      "in", " ", 
      RowBox[{"h_i", "'"}], "s"}]}]}], ",", 
   RowBox[{"where", " ", "F", 
    RowBox[{"(", "t", ")"}], " ", "is", " ", "the", " ", "instantaneous", " ",
     
    RowBox[{"frequency", "."}]}]}]}]], "Text",
 CellChangeTimes->{
  3.472933485856436*^9, {3.472934739593316*^9, 3.472934824292151*^9}, {
   3.472934996709121*^9, 3.472935025035933*^9}, {3.47293520341917*^9, 
   3.472935204338819*^9}, {3.473088563881653*^9, 3.473088588810369*^9}, {
   3.473110439776879*^9, 3.47311050681657*^9}, {3.473130648090571*^9, 
   3.473130664604451*^9}, {3.473132309310089*^9, 3.473132330894548*^9}, {
   3.47362070301891*^9, 3.473620704216636*^9}},
 FontFamily->"Bitstream Vera Sans",
 FontSize->18],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"c", " ", "=", " ", 
   RowBox[{"Cos", "[", "iota", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   SubscriptBox["a", "1"], " ", "=", " ", 
   RowBox[{
    FractionBox["1", "r"], 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"-", 
        RowBox[{"(", 
         RowBox[{"1", "+", 
          RowBox[{"c", "^", "2"}]}], ")"}]}], " ", 
       RowBox[{"Cos", "[", 
        SubscriptBox["\[CapitalPhi]", "0"], "]"}], " ", 
       RowBox[{"Sin", "[", 
        RowBox[{"2", " ", "\[Psi]"}], "]"}]}], "-", 
      RowBox[{"2", " ", "c", " ", 
       RowBox[{"Cos", "[", 
        RowBox[{"2", " ", "\[Psi]"}], "]"}], " ", 
       RowBox[{"Sin", "[", 
        SubscriptBox["\[CapitalPhi]", "0"], "]"}]}]}], ")"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   SubscriptBox["a", "2"], " ", "=", 
   RowBox[{
    FractionBox["1", "r"], 
    RowBox[{"(", " ", 
     RowBox[{
      RowBox[{"2", " ", "c", " ", 
       RowBox[{"Cos", "[", 
        SubscriptBox["\[CapitalPhi]", "0"], "]"}], " ", 
       RowBox[{"Sin", "[", 
        RowBox[{"2", " ", "\[Psi]"}], "]"}]}], "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "+", 
         RowBox[{"c", "^", "2"}]}], ")"}], " ", 
       RowBox[{"Cos", "[", 
        RowBox[{"2", " ", "\[Psi]"}], "]"}], " ", 
       RowBox[{"Sin", "[", 
        SubscriptBox["\[CapitalPhi]", "0"], "]"}]}]}], ")"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   SubscriptBox["a", "3"], " ", "=", " ", 
   RowBox[{
    FractionBox["1", "r"], 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"-", 
        RowBox[{"(", 
         RowBox[{"1", "+", 
          RowBox[{"c", "^", "2"}]}], ")"}]}], " ", 
       RowBox[{"Cos", "[", 
        RowBox[{"2", " ", "\[Psi]"}], "]"}], " ", 
       RowBox[{"Cos", "[", 
        SubscriptBox["\[CapitalPhi]", "0"], "]"}]}], "+", 
      RowBox[{"2", " ", "c", " ", 
       RowBox[{"Sin", "[", 
        RowBox[{"2", " ", "\[Psi]"}], "]"}], " ", 
       RowBox[{"Sin", "[", 
        SubscriptBox["\[CapitalPhi]", "0"], "]"}]}]}], ")"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    SubscriptBox["a", "4"], " ", "=", " ", 
    RowBox[{
     FractionBox["1", "r"], 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"-", "2"}], " ", "c", " ", 
        RowBox[{"Cos", "[", 
         RowBox[{"2", " ", "\[Psi]"}], "]"}], " ", 
        RowBox[{"Cos", "[", 
         SubscriptBox["\[CapitalPhi]", "0"], "]"}]}], "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"1", "+", 
          RowBox[{"c", "^", "2"}]}], ")"}], " ", 
        RowBox[{"Sin", "[", 
         RowBox[{"2", " ", "\[Psi]"}], "]"}], " ", 
        RowBox[{"Sin", "[", 
         SubscriptBox["\[CapitalPhi]", "0"], "]"}]}]}], ")"}]}]}], ";"}], 
  "\n"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"h1", " ", "=", " ", 
   RowBox[{"B", " ", 
    RowBox[{"(", 
     RowBox[{"2", "M", " ", "\[Eta]"}], ")"}], "x", " ", 
    RowBox[{"Cos", "[", "\[CapitalPhi]", "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"h2", " ", "=", " ", 
   RowBox[{"A", " ", 
    RowBox[{"(", 
     RowBox[{"2", "M", " ", "\[Eta]"}], ")"}], "x", " ", 
    RowBox[{"Sin", "[", "\[CapitalPhi]", "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"h3", " ", "=", " ", 
   RowBox[{"A", " ", 
    RowBox[{"(", 
     RowBox[{"2", "M", " ", "\[Eta]"}], ")"}], "x", " ", 
    RowBox[{"Cos", "[", "\[CapitalPhi]", "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"h4", " ", "=", " ", 
   RowBox[{"B", " ", 
    RowBox[{"(", 
     RowBox[{"2", "M", " ", "\[Eta]"}], ")"}], "x", " ", 
    RowBox[{"Sin", "[", "\[CapitalPhi]", "]"}]}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.473087955724536*^9, 3.473088149987*^9}, {
   3.473088211324668*^9, 3.473088237405479*^9}, {3.473088416008078*^9, 
   3.473088422085199*^9}, {3.473088532260028*^9, 3.47308855836538*^9}, {
   3.473088595030854*^9, 3.473088640827654*^9}, {3.473088950576395*^9, 
   3.47308895507785*^9}, {3.473110261814255*^9, 3.473110433549366*^9}, {
   3.473132340065735*^9, 3.473132515527796*^9}, {3.473620708822773*^9, 
   3.473620720034397*^9}, {3.473630794613918*^9, 3.473630799982151*^9}, {
   3.473630834568567*^9, 3.473630925194372*^9}, {3.473632782038197*^9, 
   3.473632786613661*^9}, 3.474043618137027*^9, {3.474043800856585*^9, 
   3.47404380354937*^9}, {3.47405410134217*^9, 3.474054118939833*^9}},
 FontSize->16],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "-", 
      SuperscriptBox[
       RowBox[{"Cos", "[", "iota", "]"}], "2"]}], ")"}], " ", 
    RowBox[{"Cos", "[", 
     SubscriptBox["\[CapitalPhi]", "0"], "]"}], " ", 
    RowBox[{"Sin", "[", 
     RowBox[{"2", " ", "\[Psi]"}], "]"}]}], "-", 
   RowBox[{"2", " ", 
    RowBox[{"Cos", "[", "iota", "]"}], " ", 
    RowBox[{"Cos", "[", 
     RowBox[{"2", " ", "\[Psi]"}], "]"}], " ", 
    RowBox[{"Sin", "[", 
     SubscriptBox["\[CapitalPhi]", "0"], "]"}]}]}], "r"]], "Output",
 CellChangeTimes->{3.474054112186159*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"I1", " ", "=", " ", 
   RowBox[{
    SubsuperscriptBox["\[Integral]", "0", "Pi"], " ", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         SubscriptBox["a", "1"], " ", "B"}], " ", "+", " ", 
        RowBox[{
         SubscriptBox["a", "3"], " ", "A"}]}], ")"}], "^", "2"}], " ", 
     RowBox[{"\[DifferentialD]", "\[Psi]"}]}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"I2", " ", "=", " ", 
   RowBox[{
    SubsuperscriptBox["\[Integral]", 
     RowBox[{"-", "Pi"}], "Pi"], 
    RowBox[{"I1", 
     RowBox[{"\[DifferentialD]", 
      SubscriptBox["\[CapitalPhi]", "0"]}]}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Simplify", "[", 
  RowBox[{"I2", "/", 
   RowBox[{"(", 
    RowBox[{"2", " ", 
     RowBox[{"Pi", "^", "2"}]}], ")"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.4740542271525*^9, 3.47405442744947*^9}, {
  3.474054589543489*^9, 3.474054592902991*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     SuperscriptBox["A", "2"], "+", 
     SuperscriptBox["B", "2"]}], ")"}], " ", 
   RowBox[{"(", 
    RowBox[{"35", "+", 
     RowBox[{"28", " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "iota"}], "]"}]}], "+", 
     RowBox[{"Cos", "[", 
      RowBox[{"4", " ", "iota"}], "]"}]}], ")"}]}], 
  RowBox[{"32", " ", 
   SuperscriptBox["r", "2"]}]]], "Output",
 CellChangeTimes->{3.474054259174334*^9, 3.474054345195495*^9, 
  3.47405440035414*^9, 3.474054440002543*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"I1", " ", "=", " ", 
   RowBox[{
    SubsuperscriptBox["\[Integral]", "0", 
     RowBox[{"Pi", "/", "2"}]], 
    RowBox[{
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         SubscriptBox["a", "1"], "B"}], "  ", "+", " ", 
        RowBox[{
         SubscriptBox["a", "3"], "A"}]}], ")"}], "2"], 
     RowBox[{"\[DifferentialD]", " ", "iota"}]}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"I2", " ", "=", " ", 
   RowBox[{
    SubsuperscriptBox["\[Integral]", "0", "Pi"], 
    RowBox[{"I1", 
     RowBox[{"\[DifferentialD]", "\[Psi]"}]}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"I3", " ", "=", " ", 
   RowBox[{
    SubsuperscriptBox["\[Integral]", 
     RowBox[{"-", "Pi"}], "Pi"], 
    RowBox[{"I2", 
     RowBox[{"\[DifferentialD]", 
      SubscriptBox["\[CapitalPhi]", "0"]}]}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Simplify", "[", 
  RowBox[{"I3", "/", 
   RowBox[{"(", 
    RowBox[{"Pi", "^", "3"}], ")"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.473619941000845*^9, 3.473620092678078*^9}, {
   3.473620151531045*^9, 3.473620415774559*^9}, {3.473620462131972*^9, 
   3.473620489402581*^9}, {3.473620525700101*^9, 3.473620562031645*^9}, {
   3.473620724919126*^9, 3.473620725094134*^9}, {3.473620878368594*^9, 
   3.473620917406889*^9}, {3.473620987187619*^9, 3.473621063450031*^9}, {
   3.473623758392454*^9, 3.473623792710865*^9}, {3.473623857084116*^9, 
   3.473623927724914*^9}, {3.473623958659821*^9, 3.473623961907007*^9}, {
   3.473624238559912*^9, 3.473624313320732*^9}, {3.473630036337952*^9, 
   3.473630045255339*^9}, 3.473630086321511*^9, {3.473630322542761*^9, 
   3.473630333013413*^9}, 3.473630484783049*^9, {3.473630913348974*^9, 
   3.473630913895568*^9}, {3.473630968475529*^9, 3.473631119987483*^9}, {
   3.473631184319066*^9, 3.473631236810354*^9}, 3.473631281554985*^9, {
   3.473631327727201*^9, 3.473631348399885*^9}, {3.473632812240335*^9, 
   3.473632824331852*^9}, {3.473632866937539*^9, 3.473632899883797*^9}, {
   3.473632934755433*^9, 3.473632935278668*^9}, {3.473634118035565*^9, 
   3.473634122656835*^9}, {3.474053773976148*^9, 3.474053822226861*^9}, {
   3.474053853442127*^9, 3.474053872139718*^9}, {3.474053903502026*^9, 
   3.474053951998161*^9}, {3.474054009717188*^9, 3.474054016491828*^9}, {
   3.474054126321851*^9, 3.474054145191127*^9}, 3.47405421012967*^9, {
   3.474054508795624*^9, 3.474054511750331*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{"35", " ", 
   RowBox[{"(", 
    RowBox[{
     SuperscriptBox["A", "2"], "+", 
     SuperscriptBox["B", "2"]}], ")"}]}], 
  RowBox[{"32", " ", 
   SuperscriptBox["r", "2"]}]]], "Output",
 CellChangeTimes->{
  3.473631065374536*^9, 3.473631121299918*^9, {3.473631188639843*^9, 
   3.473631237883002*^9}, 3.473631282715692*^9, 3.473631350841691*^9, 
   3.473632760686508*^9, {3.473632807094459*^9, 3.473632814190037*^9}, 
   3.473632845779782*^9, {3.473632888895123*^9, 3.473632900999704*^9}, 
   3.473632941317494*^9, 3.473634128039911*^9, {3.474053823347665*^9, 
   3.474053831855096*^9}, 3.47405387291705*^9, 3.474053904717307*^9, 
   3.474053953835006*^9, 3.474054176494499*^9, 3.474054522150867*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"I1", " ", "=", " ", 
   RowBox[{
    SubsuperscriptBox["\[Integral]", "0", "Pi"], " ", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         SubscriptBox["a", "4"], " ", "B"}], " ", "+", " ", 
        RowBox[{
         SubscriptBox["a", "2"], " ", "A"}]}], ")"}], "^", "2"}], " ", 
     RowBox[{"\[DifferentialD]", "\[Psi]"}]}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"I2", " ", "=", " ", 
   RowBox[{
    SubsuperscriptBox["\[Integral]", 
     RowBox[{"-", "Pi"}], "Pi"], 
    RowBox[{"I1", 
     RowBox[{"\[DifferentialD]", 
      SubscriptBox["\[CapitalPhi]", "0"]}]}]}]}], ";"}], "\n", 
 RowBox[{"Simplify", "[", 
  RowBox[{"I2", "/", 
   RowBox[{"(", 
    RowBox[{"2", " ", 
     RowBox[{"Pi", "^", "2"}]}], ")"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.47405453098783*^9, 3.47405455663353*^9}, {
  3.474054599588216*^9, 3.474054602415342*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     SuperscriptBox["A", "2"], "+", 
     SuperscriptBox["B", "2"]}], ")"}], " ", 
   RowBox[{"(", 
    RowBox[{"35", "+", 
     RowBox[{"28", " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "iota"}], "]"}]}], "+", 
     RowBox[{"Cos", "[", 
      RowBox[{"4", " ", "iota"}], "]"}]}], ")"}]}], 
  RowBox[{"32", " ", 
   SuperscriptBox["r", "2"]}]]], "Output",
 CellChangeTimes->{3.474054617218382*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"I1", " ", "=", " ", 
   RowBox[{
    SubsuperscriptBox["\[Integral]", "0", 
     RowBox[{"Pi", "/", "2"}]], " ", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         SubscriptBox["a", "4"], " ", "B"}], " ", "+", " ", 
        RowBox[{
         SubscriptBox["a", "2"], " ", "A"}]}], ")"}], "^", "2"}], " ", 
     RowBox[{"\[DifferentialD]", "iota"}]}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"I2", " ", "=", " ", 
   RowBox[{
    SubsuperscriptBox["\[Integral]", "0", "Pi"], 
    RowBox[{"I1", 
     RowBox[{"\[DifferentialD]", "\[Psi]"}]}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"I3", " ", "=", " ", 
   RowBox[{
    SubsuperscriptBox["\[Integral]", 
     RowBox[{"-", "Pi"}], "Pi"], 
    RowBox[{"I2", 
     RowBox[{"\[DifferentialD]", 
      SubscriptBox["\[CapitalPhi]", "0"]}]}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"Simplify", "[", 
   RowBox[{"I3", "/", 
    RowBox[{"(", 
     RowBox[{"Pi", "^", "3"}], ")"}]}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.474054658122461*^9, 3.474054710719403*^9}, {
  3.474059792981289*^9, 3.474059793141175*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{"35", " ", 
   RowBox[{"(", 
    RowBox[{
     SuperscriptBox["A", "2"], "+", 
     SuperscriptBox["B", "2"]}], ")"}]}], 
  RowBox[{"32", " ", 
   SuperscriptBox["r", "2"]}]]], "Output",
 CellChangeTimes->{3.474054723535375*^9, 3.474059806207314*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"P", " ", "=", " ", 
   RowBox[{
    SqrtBox[
     FractionBox[
      RowBox[{"5", " ", "eta"}], "96"]], 
    SuperscriptBox["M", 
     RowBox[{"5", "/", "6"}]], 
    SuperscriptBox["\[Pi]", 
     RowBox[{
      RowBox[{"-", "2"}], "/", "3"}]], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"a1", " ", "B"}], " ", "+", " ", 
      RowBox[{"a3", " ", "A"}]}], ")"}]}]}], " ", 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Q", " ", "=", " ", 
    RowBox[{
     RowBox[{"-", 
      SqrtBox[
       FractionBox[
        RowBox[{"5", " ", "eta"}], "96"]]}], 
     SuperscriptBox["M", 
      RowBox[{"5", "/", "6"}]], 
     SuperscriptBox["\[Pi]", 
      RowBox[{
       RowBox[{"-", "2"}], "/", "3"}]], 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"a4", " ", "B"}], " ", "+", " ", 
       RowBox[{"a2", " ", "A"}]}], ")"}]}]}], " ", ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{"FullSimplify", "[", 
  RowBox[{
   RowBox[{"D", "[", 
    RowBox[{"Q", ",", "a4"}], "]"}], "/", "Q"}], "]"}]}], "Input",
 CellChangeTimes->{{3.474063476995604*^9, 3.474063478412742*^9}, {
  3.474063599735939*^9, 3.474063603389741*^9}, {3.474063656854558*^9, 
  3.47406400464851*^9}, {3.474064059221933*^9, 3.474064072837292*^9}, {
  3.474231589791453*^9, 3.474231594983514*^9}, {3.474399465407821*^9, 
  3.474399501783907*^9}, {3.474400550953182*^9, 3.474400552033049*^9}, {
  3.474400647175838*^9, 3.474400647789782*^9}, {3.474400764411647*^9, 
  3.474400765164373*^9}, {3.474400953689777*^9, 3.474400954422452*^9}, {
  3.474401025323604*^9, 3.474401025826267*^9}, {3.4744010938544*^9, 
  3.474401116711491*^9}, {3.474401454906023*^9, 3.47440146277887*^9}, {
  3.47440156547543*^9, 3.47440158919314*^9}, {3.474401629604264*^9, 
  3.474401630162782*^9}, {3.474401696447822*^9, 3.474401700062605*^9}, {
  3.474401740561548*^9, 3.474401741378049*^9}}],

Cell[BoxData[
 FractionBox["B", 
  RowBox[{
   RowBox[{"A", " ", "a2"}], "+", 
   RowBox[{"a4", " ", "B"}]}]]], "Output",
 CellChangeTimes->{{3.474399473479786*^9, 3.474399502339897*^9}, 
   3.474400552885462*^9, 3.474400654977944*^9, 3.474400765701967*^9, 
   3.474400955115709*^9, 3.474401027177155*^9, {3.47440109497873*^9, 
   3.474401117351952*^9}, {3.474401456109023*^9, 3.47440146327094*^9}, {
   3.47440157482423*^9, 3.474401589677763*^9}, 3.474401630767658*^9, {
   3.474401697568698*^9, 3.47440170051263*^9}, 3.474401741861579*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"PAngAvg", " ", "=", " ", 
     RowBox[{
      SqrtBox[
       FractionBox[
        RowBox[{"35", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"A", "^", "2"}], " ", "+", " ", 
           RowBox[{"B", "^", "2"}]}], ")"}]}], "32"]], 
      SqrtBox[
       FractionBox[
        RowBox[{"5", " ", "eta"}], "96"]], 
      FractionBox[
       RowBox[{
        RowBox[{"M", "^", 
         RowBox[{"(", 
          RowBox[{"5", "/", "6"}], ")"}]}], " ", 
        SuperscriptBox["\[Pi]", 
         RowBox[{
          RowBox[{"-", "2"}], "/", "3"}]]}], "r"]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"FullSimplify", "[", "PAngAvg", "]"}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"FullSimplify", "[", 
    RowBox[{
     RowBox[{"D", "[", 
      RowBox[{
       RowBox[{"-", "PAngAvg"}], ",", " ", "B"}], "]"}], "/", 
     RowBox[{"-", "PAngAvg"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.47439868302299*^9, 3.474398779636456*^9}, {
  3.47439883385472*^9, 3.474398834998109*^9}, {3.474398874283882*^9, 
  3.474398875194704*^9}, {3.474398927085002*^9, 3.474398929502382*^9}, {
  3.474398963022785*^9, 3.474398963248197*^9}, {3.474399010347177*^9, 
  3.474399010715105*^9}, {3.474399263606957*^9, 3.47439927531342*^9}, {
  3.474399326302479*^9, 3.474399357799383*^9}, {3.474399410499304*^9, 
  3.474399442659846*^9}, {3.474400542623446*^9, 3.474400543078748*^9}, {
  3.474400658598158*^9, 3.474400659781198*^9}, {3.474400768972881*^9, 
  3.474400769619853*^9}, {3.474400976414973*^9, 3.474401000078864*^9}, {
  3.47440147551315*^9, 3.474401483874853*^9}, {3.474401578568013*^9, 
  3.474401593706024*^9}, {3.474401625179455*^9, 3.474401625570744*^9}}],

Cell[BoxData[
 FractionBox["B", 
  RowBox[{
   SuperscriptBox["A", "2"], "+", 
   SuperscriptBox["B", "2"]}]]], "Output",
 CellChangeTimes->{{3.474398767438095*^9, 3.474398780127885*^9}, 
   3.474398835821449*^9, 3.474398875926698*^9, 3.474398930559196*^9, 
   3.474398963975875*^9, 3.474399011292066*^9, 3.474399275827468*^9, {
   3.47439932815955*^9, 3.474399358372908*^9}, {3.47439941207374*^9, 
   3.4743994439555*^9}, 3.474400548252905*^9, 3.474400660300291*^9, 
   3.474400770242277*^9, {3.47440098136088*^9, 3.47440100043933*^9}, {
   3.474401480751452*^9, 3.474401484324644*^9}, {3.474401579682401*^9, 
   3.474401594298426*^9}, 3.474401626043385*^9}]
}, Open  ]]
},
WindowSize->{1258, 875},
WindowMargins->{{41, Automatic}, {Automatic, 53}},
ShowSelection->True,
FrontEndVersion->"6.0 for Linux x86 (32-bit) (June 19, 2007)",
StyleDefinitions->FrontEnd`FileName[{"Creative"}, "NaturalColor.nb", 
  CharacterEncoding -> "iso8859-1"]
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[568, 21, 1060, 17, 65, "Text"],
Cell[1631, 40, 1595, 36, 74, "Input"],
Cell[CellGroupData[{
Cell[3251, 80, 553, 16, 65, "Input"],
Cell[3807, 98, 657, 21, 52, "Output"]
}, Open  ]],
Cell[4479, 122, 1080, 17, 65, "Text"],
Cell[5562, 141, 1863, 40, 74, "Input"],
Cell[7428, 183, 1040, 17, 65, "Text"],
Cell[CellGroupData[{
Cell[8493, 204, 1746, 37, 46, "Input"],
Cell[10242, 243, 2705, 78, 76, "Output"]
}, Open  ]],
Cell[12962, 324, 1633, 42, 109, "Text"],
Cell[CellGroupData[{
Cell[14620, 370, 4482, 129, 418, "Input"],
Cell[19105, 501, 629, 19, 61, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[19771, 525, 962, 30, 119, "Input"],
Cell[20736, 557, 548, 17, 62, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[21321, 579, 2468, 55, 157, "Input"],
Cell[23792, 636, 740, 16, 62, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[24569, 657, 943, 29, 119, "Input"],
Cell[25515, 688, 480, 16, 62, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26032, 709, 1127, 35, 157, "Input"],
Cell[27162, 746, 288, 9, 62, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[27487, 760, 1913, 50, 193, "Input"],
Cell[29403, 812, 541, 10, 55, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[29981, 827, 1767, 44, 152, "Input"],
Cell[31751, 873, 659, 12, 56, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
